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ABSTRACT 


An accreting black hole is, by definition, characterized by the drain. Namely, 
the matter falls into a black hole much the same way as water disappears down 
a drain - matter goes in and nothing comes out. As this can only happen in 
a black hole, it provides a way to see “a black hole”, an unique observational 
signature of black holes. The accretion proceeds almost in a free-fall manner 
close to the black hole horizon, where the strong gravitational field dominates the 
pressure forces. In this paper we calculate (by using Monte-Carlo simulations) 
the specific features of X-ray spectra formed as a result of upscattering of the 
soft (disk) photons in the converging inflow (Cl) within about 3 Schwarzschild 
radii of the black hole. The full relativistic treatment has been implemented to 
reproduce these spectra. We show that spectra in the soft state of black hole 
systems (BHS) can be described as the sum of a thermal (disk) component and 
the convolution of some fraction of this component with the Cl upscattering 
spread (Green’s) function. The latter boosted photon component is seen as 
an extended power-law at energies much higher than the characteristic energy 
of the soft photons. We demonstrate the stability of the power spectral index 
{a = 1.8 ±0.1) over a wide range of the plasma temperature 0 — 10 keV and mass 
accretion rates (higher than 2 in Eddington units). We also demonstrate that 
the sharp high energy cutoff occurs at energies of 200-400 keV which are related 
to the average energy of electrons rrieC^ impinging upon the event horizon. 

The spectrum is practically identical to the standard thermal Comptonization 
spectrum (Hua & Titarchuk 1995) when the Cl plasma temperature is getting 
of order of 50 keV (the typical ones for the hard state of BHS). In this case one 
can see the effect of the bulk motion only at high energies where there is an 
excess in the Cl spectrum with respect to the pure thermal one. Furthermore 
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we demonstrate that the change of spectral shapes from the soft X-ray state 
to the hard X-ray state is clearly to be related with the temperature of the 
bulk flow. In other words the effect of the bulk Comptonization compared to 
the thermal one is getting stronger when the plasma temperature drops below 
10 keV. These Monte-Carlo simulated Cl spectra are a inevitable stamp of the 
BHS where the strong gravitational held dominates the pressure forces. 

Subject headings: Black hole physics - X-rays: general - radiation 
mechanism: non thermal 
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1. Introduction 


A lot of theories and models have been suggested to explain the observational situations 
of black hole systems (BHS); for recent reviews of the astrophysics of black holes [see e.g., 
Liang fc Narayan (1997)| , [Zhang et al. (1997]|, [Titarchuk (1996)|, [Titarchuk (1997) . 


The BHS are generally observed in two spectral states: one with a high luminosity 
soft thermal bump and a low luminosity power law extending to several hundred of keV 
(hereafter called the “soft state”), and the other at high luminosity in high energies (> 10 
keV) showing a typical thermal Comptonization component (called the “hard state” all 
along the paper). 


It is now clear that the observed X-ray spectra resulting from the accretion of matter 
onto the central object is divided according to the connection between the two effects: (i). 
The gravitational energy of matter is released in the disk due to the viscous dissipation 
between the adjacent layers (Shakura & Sunyaev 1973, hereafter SS73). The energy release 
occurs in the optically thick medium and thus the disk emits a soft blackbody like radiation 
(related to the effective disk area, SS73) of a fraction of keV and tens of eV for galactic and 
extragalatic sources respectively, (ii). Another mechanism is related with the geometric 
compression of the matter accreting in advection dominated manner where the radial 
velocity is some fraction of the free fall velocity but the viscous heating can be essential 
in the flow too (e.g. Chakrabarti & Titarchuk 1995, hereafter CT95; Narayan & Yi 1994). 
Energy is mostly transferred by the protons due to their much higher mass. These protons 
afterward transfer their energy to the ambient electrons through Coulomb collisions in a 
relatively dense cloud (called the Compton cloud) close to the black hole. This cloud can 
be a result of the propagation (accretion) of the sub-Keplerian component towards the 
black hole. In such a type of accretion, shocks or centrifugal barrier regions can be formed 
(CT95). 
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Since the cross-sections for radiative processes are inversely proportional to the mass 
of the emitter, it is the electrons which then radiatively cool the plasma. In such an 
high temperature plasma, the main electron cooling process is the Comptonization on the 
ambient soft photons. The plasma temperature of the zone emitting the hard radiation is 
regulated by the supply of soft photons from the disk. 

It was demonstrated in CT95 that the hard-soft state transition in BHS is regulated by 
redistribution of mass accretion rates between the quasi-spherical (sub-Keplerian) and the 
Keplerian disk components. It was also shown there for the first time that the converging 
inflow spectral feature (as an intrinsic signature of black hole) can be seen in the soft 
state. Namely, near the horizon, the strong gravitational held is expected to dominate the 
pressure forces and thus to drive the accreting material in almost a free fall i.e. a converging 
inhow (Titarchuk & Zannias 1998, hereafter TZ98). Unlike, for other compact objects the 
pressure forces are becoming dominant as their surface is approached, and thus such a free 
falling state is absent. In contrast CT95 emphasized that in the soft state, the optically 
thick Compton cloud surrounding a black hole is cooled down by the disk soft photons and 
the extended power law will be observed. 

The fact of the presence of the geometrically thin disk in BH system leads to the very 
important statement regarding the dominance of gravitational forces over radiative forces: 
Because the gravitational force in the disk is a small H/R fraction of the gravitational force 
acting in the radial direction (where H and R are the disk height and radius respectively) the 
local luminosity in the disk cannot be higher than H/R times the Eddington limit. At some 
certain point the disk is disrupted or terminated, for example due to General Relativistic 
effects: the stable orbits do not exist below three Schwarzschild radii. Then the matter is 
forced to proceed towards the black hole in almost free-fall manner - the radiative forces 
being negligible with respect to the gravitational forces. The next question now is to know 
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how the soft disk photons illuminating the bulk inflow can be upcomptonized to energies of 
the order of mgC^. 

Titarchuk & Zannias 1998 show that the production of hard photons as a function of 
the radius from the black hole peaks at around 2 Schwarzschild radii {2rs) and most of the 
hard continuum photons observed are produced within Sr* from the black hole. Shrader & 
Titarchuk 1998 (hereafter ShT98) demonstrate that the emergent spectra of the converging 
inflow atmosphere successfully explain the continuum X-ray spectra of black hole systems 
(BHS) in their soft state. 

To check these statements and predictions we have simulated the emergent energy 
spectrum, treating as precisely as possible the Comptonization process. To do this, we need 
to know hrst what is the movement of the high energy electrons: it is in fact, a composition 
of the free-fall motion onto the black hole and of the Brownian (thermal) motion which 
can be characterized by a temperature. Whereas the thermal Comptonization has been 
already explored in several ways (e.g. Sunyaev & Titarchuk 1980; Pozdnyakov, Sobol’ & 
Sunyaev 1983, hereafter PSS; Nagirner & Poutanen 1994; Titarchuk 1994, hereafter T94; 
Titarchuk & Lyubarskij 1995; Hua & Titarchuk 1995; Giesler & Kirk 1997), the effects of 
the free-fall motion of the electrons on the resulting spectrum is not yet clearly understood. 
Also, the general relativistic effects close to the black hole were not yet plainly taken into 
account. In series of papers, Titarchuk, Mastichiadis & Kylafis (1996), (1997) (hereafter 
TMK96 and TMK97) present the exact numerical and approximate analytical solutions 
of the problem of spectral formation in a converging inflow, taking into account the inner 
boundary condition, the dynamical effects of the free fall, and the thermal motion of the 
electrons. The inner boundary has been taken at a finite radius with the spherical surface 
being considered as a fully absorptive. 

TMK have used a variant of the Fokker-Plank formalism where the inner boundary 
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mimics a black-hole horizon; no relativistic effects (special or general) have been taken into 
account. Thus their results are instructively useful but they are not directly comparable with 
the observations. But using these numerieal and analytical techniques they demonstrated 
that the extended power laws are present in the resulting spectra in addition to the blackbody 
like emission at lower energies. 

Our approach is to determine a realistic emergent Comptonization spectrum using 
two different approaches: one, purely analytical, is described in detail in TZ98, and 
the other, numerical, is presented below. TZ98 analyzed the exact general relativistic 
integro-differential equation of radiative transfer describing the interaction of low energy 
photons with a Maxwellian distribution of the relativistic electrons in the gravitational 
held of a Schwarzschild black hole. They proved that due to Comptonization an initial 
arbitrary spectrum of low energy photons unavoidably results in spectra characterized by an 
extended power-law feature. They examined the spectral index a by using both analytical 
and numerical methods in the lower temperature limit (0 = kTe/meC^ -C 1), where Te is 
the electron temperature and meC^ the electron rest mass energy, and with a varying mass 
accretion rate. They demonstrated the stability of the {a ~ 1.8) power-law, which is mainly 
due to the asymptotic independence of the spectral index on the mass accretion rate and 
its weak dependence on plasma temperatures in the lower temperature limit. In fact, all 
corrections to the index due the temperature effects are related with the corrections of the 
kernel of the relativistic kinetic equation (TZ98) which are of the order of the dimensionless 
temperature 0. 

In this paper, we will describe the numerical approach briefly outlining the main 
features of the Monte-Carlo simulation code in §2. We will present the results in the “Flat” 
case in §3 for which we assume the propagation of the photons in straight lines. Then, we 
will show what we have obtained in the most realistic general relativistic case in §4. In 



sections 3-4, we will compare also these results to the ones obtained by different (analytical 
and numerical) methods in the same physical framework (TMK97 and TZ98). We will 
show that the observed spectral transitions in BHS is clearly to be related to the bulk flow 
temperature in §5. We will demonstrate the relevance of the converging inflow model to 
the recent high energy observations in §6. Finally, we will summaries our work and draw 
conclusions in 57. 


2. The Monte-Carlo simulation 


As explained above, this signature originates from upscattering of low energy photons 
by fast moving electrons with velocities, n, approaching the speed of light, c. A soft photon 
of energy E, in the process of multiple scattering off the electrons, gets substantially 
blue-shifted. For a single scattering, the photon scattered off an electron moving with 
velocity v, changes its energy due to Doppler effect as follows 

1 — (n/c) cos 6* 


E' = E- 


( 1 ) 


1 — {v/c) cos 9' 

The photon is then significantly blue-shifted if, at least, it gets once scattered in the 
direction of the electron motion (i.e. when cos 6' ~ 1). This is schematically illustrated 
in Figure 1 for two scattering events. In the first event the soft photon is effectively 
upcomptonized and in the second one the photon is scattered in such a direction that it can 
be detected by an observer. 


Editor, please put Figure 1 here. 

In the first scattering event we assume the direction of incident photon, 6i, to be 
nearly normal to the electron velocity, and the direction of the scattered photon to be 
nearly aligned with the electron velocity. During its outward propagation through the 
converging-inflow medium, the angle between the photon and electron velocity increases 
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Thus, in the second event the cosine angle, cos 02, tends to approach zero. The angle of 
outgoing photon, 62 , has to be large enough, in order for the Doppler boosted photon to 
reach an observer. 


The geometry mostly used in these simulations is shown in Figure 2: it consists of a thin 
disk surronnding a spherical cloud harbouring a black hole in its center. The disk is always 
supposed to be optically thick, whereas we are assuming a free fall for the backgronnd flow 
where the bulk velocity of the infalling plasma is given by v{r) = c{rs/ry^‘^- 

Editor, please put Figure 2 here. 


In our simulations (see also the kinetic formalism in TZ98) we use the number 
density n measured in the local rest frame of the flow that is n = m{rs/rY^‘^/{ 2 raT)- 
Here m = M/Me, M is the mass accretion rate, aT is the Thomson cross section. 

Me = Le/c^ = 4:7rGMmp/aTC is the Eddington accretion rate and M is the black hole 
mass. Then the optical thickness of the clond is compnted through the mass accretion rate 
771 according to the following formulae : 



where Vgut and r* are the onter boundary and the Schwarzschild radii, respectively. In most 
of all results given below, except otherwise presented in part 4, the clond radius rout has 
been taken to be 3 r*. The different values of m which has been used in those simulations 


are tti = 0.5,1,1.5, 2,4, 7. 


We have varied the clond electron temperatnre from 0 to 50 keV. In addition to 
the thermal motion of the electrons in the cloud, we have also taken into acconnt their 
free-fall onto the central black hole. The seed X-ray photons were generated uniformly and 
isotropically at the snrface of the inner edge of the accretion border of the disk. For our 
simulations we use as an example the disk thermal spectrum with a temperature of 0.5 keV. 
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In the simulations, we have followed the trajectory of each photons in the following 
way: hrst, we have made a random generation according to an exponential law: P{t) = e~'^ 
to get the value of the optical thickness T that a photon will travel before interacting. 
Then we have integrated the optical thickness r = pal along the photon path up to T, 
taking into account the variation of the cloud density p and of the cross section a with the 
radius. The gravitational shift endured by the photon was also computed at each step of 
this integration. 

The photon paths are straight lines in what we called the “Flat” case, or curves in 
terms of the Schwarzschild geometry in the general relativistic case. In fact, the variation of 
the photon energy, and of the photon position (with respect to the local normal) along the 
photon trajectory can be extracted from the full kinetic equation presented in TZ98 (Eq. 
12). The characteristics of the differential part of the equation provides these dependencies. 

If, at the end of the integration, the photon has not left the cloud, we simulate a 
Compton scattering with an electron according to the method described in PSS, taking 
into account the exact motion of the electron, that is a composition of its free-fall motion 
onto the black hole with its Brownian thermal motion. To do this, we hrst compute the 
scattering electron momentum, and derive the scattered photon and electron characteristics 
from the Compton scattering kinetics. We then check if this event is consistent with the 
Compton scattering probabilities; if yes the event is kept; if not another scattering electron 
is generated, and the process goes on until the event is accepted. This process has been 
successfully checked by comparing its results with the analytical ones of Hua & Titarchuk 
1995 in the case of relativistic thermal electrons. 

Once the new energy and direction of the photon been determined by the Compton 
kinematics, we track it in the same way as above until it makes another scattering, or it 
escapes from the cloud, or until it is “absorbed” by the black hole at its horizon. 
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The parameters of the simulations are the Compton cloud electron temperature, 
the cloud outer radius, Tout and the mass accretion rate, rh. 


3. “Flat” Case 

We have hrst investigated what could be the emergent spectrum in the “Flat” case, 
that is supposing that photons goes in straight lines. But the relativistic effects such 
as gravitational redshift, dependence of the Klein-Nishina cross section on energy, or 
“absorption” of photons at the black hole horizon were correctly taken into account in these 
computations. Even if this case is not physically consistent, as we cannot separate in theory 
the gravitational redshift effects from the light bending, it is instructive as it will enable 
us to derive the consequences of space curvature around the black hole on the emerging 
spectrum, by comparing this case to the general relativistic one described below. 

We have calculated the spectrum emerging from a plasma cloud of temperature 5 keV 
for rh equal to 2 with an outer radius. Tout = Sr^. They are typically the physical conditions 
of the soft state. As foreseen theoretically, (TMK96-97) the spectrum observed at inhnity 
consist of a soft component (blackbody like) coming from input photons that escaped after 
a few scattering without any signihcant energy change and the high energy part (a power 
law) extending to high energies of order 300 keV. The shape of spectrum is reminiscent 
of what is observed from BHS in their soft state (ShT98) and it is almost identical to the 
spectrum presented in TMK97 (Fig. 1). 

However it is clear that even in the low temperature case, the spectrum we obtained 
is too hard to account for what is observed from soft state BHS. We have indeed htted the 
Monte-Carlo spectra (see below Fig. 5) with a test spectrum, composed of two components: 
a Planckian spectrum with a hxed temperature of the soft photons and convolution of some 
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fraction of this component with the converging inflow upscattering Green’s function. The ht 
parameters were then the color temperature of the soft component, the normalization of the 
upscattering component relative to the soft one, and the spectral index of the power law. 
An example of this £t to the spectrum shown in Figure 5 (which is a spectrum obtained in 
the General Relativistic case, with kT^, = 20 keV and m = 4), is given in Figure 3. 

The fitting parameters are kTc = 0.4 keV, the fraction of the illumination of the bulk 
inflow atmosphere by the disk photons, / = 0.48, and the spectral index a = 1.5. It is worth 
noting that the htting temperature kT^ = 0.4 keV measured by an observer is less than the 
input temperature of the soft photons 0.5 keV because of the gravitational redshift. 

Editor, please put figure 3 here. 

So, as it could be seen in Table 1, the spectral index we found by fitting the computed 
spectrum in the “Flat” case, are around or harder than 1, which is harder than it has been 
observed in the soft state (around 1.5, see for instance Ebisawa, Titarchuk & Chakrabarti, 
1996, hereafter ETC). It is worth noting that the spectral indices presented in table 1, 
are very close to TMK97 results (see Figs. 6 there) made with an assumption of a purely 
Newtonian geometry. 

So, as expected, the modelling we made in the “Flat” case is not realistic. It has been 
shown moreover that it is possible to soften the emergent spectrum by taking into account 
general relativistic curvature effects (TZ98). We have verified this result by computing the 
Gomptonized spectrum in a fully general relativistic case, taking into account the space 
curvature around the black hole horizon. 
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4. The General Relativistic Case 

For this case we have modihed our simulation code in order to take into account the 
space curvature around the black hole horizon: instead of following the photon path along a 
straight line, we have computed at each step the exact curved trajectory of the photon. The 
other physical effects (Compton diffusion, redshift, electrons thermal and free fall motion, 
etc, ...), has been computed in the same way as described previously. This description is 
then now fully consistent with the General Relativity theory. We can see in Figure 4 an 
example of the resulting emerging spectrum for kTf, = 5 keV, and m = 2: as expected, this 
spectrum is softer than the one we obtained in the “Flat” case (compare table 1 and 2 and 
also see TMK97). 

Editor, please put fig. 4 here 

A summary of the results we calculated for different physical conditions is given in 
table 2, where we have put the spectral index of the power law, while we fit our spectrum 
with a two-component test spectrum as described above. The relative error on the quoted 
value of the spectral index in table 2 is always around ±5%. We can notice from this table 
that, as discussed in TZ98, the spectral index we found in the low kT^ regime, are only 
slightly dependent of the mass accretion rate m, provided that m is larger than 1. 

Finally, we can see in table 2 that the index of the computed spectra in the soft state 
are now between 1.5 and 2.0, which is closest to the observations. 

It is worth noting here that the observed power law index of Narrow Line Seyfert 1 
(NLSl) galaxy population - which may represent the extragalactic analogue of the BHS in 
the high-soft state - is around 1.8. This is close to what we have obtained in the low kT^ 
regime, as it is expected from Comptonization on electrons whose motion is dominated by 
the bulk free-fall. 



-14- 


Table 1: Spectral index of the power law tail, obtained in the “Flat” case. 


kTe 

rh = 0.5 

171 = 1 

rh = 1.5 

rh = 2 

rh = 4 

rh = 7 

0 keV 

2.4* 

1.6* 

0.9 

0.8 

0.9 

1.2 

5 keV 

2.1* 

1.1 

0.8 

0.8 

0.9 

1.0 

10 keV 

2.0* 

1.1 

0.8 

0.8 

0.9 

1.0 

50 keV 

1.6 

1.1 

0.8 

0.9 

0.7 

0.5 


* The spectral index determination is less accurate for these values of parameters because 
the hard component of the spectrum is less pronounced in the £t region (10-50 keV). 


Table 2: Spectral index of the power law tail, obtained in the general relativistic case. 


kTe 

rh = 0.5 

rh = 1 

rh = 1.5 

rh = 2 

rh = 4 

rh = 7 

0 keV 

2.8 

2.3 

2.1 

2.05 

2.0 

1.9 

5 keV 

2.8 

2.2 

2.0 

1.9 

1.8 

1.8 

10 keV 

2.8 

2.1 

1.9 

1.9 

1.7 

1.7 

50 keV 

1.8 

1.4 

1.2 

1.0 

0.9 

0.7 
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As is seen from Figs 4-5, the power-law part of the spectrum occurs at energies E, 
lower than that of the exponential cutoff, E^. (where E^. is the average electron energy). 

Editor, please put fig. 5 here 

TZ98 have exactly solved the radiative transfer equations in the Doppler regime (when 
the photon energy is less then the energy of the electrons) in the same fully relativistic 
framework than the one we present here. Namely, they demonstrated that the Green’s 
function is a broken power law and they calculated exactly the high energy power index in 
the low temperature limit, for 0 1. 

We have compared our results and results of TZ98 obtained for zero temperature for 3 
different geometry of the Monte-Carlo simulation. The geometry used are the followings: 

1) a sphere with a monochromatic soft X-ray source (Ep = 0.5 keV) placed at the 
edge of the accretion disk, as shown in Figure 2. 

2) The same as above but with a sphere radius of 20rs. 

3) a 20r5 Comptonization cloud, with the source of 0.5 keV X-ray photons placed inside 
in a ring at hr*. 

Firstly we found only small differences in the spectral indices for all three cases, i.e. 
differences which are all within the error range of the spectral index determination, being 
around 0.1 for m greater than 2-3. Secondly, we found the mean spectral index value is 
exactly 1.8 in agreement with TZ98. We also note from this comparison that, as it is 
theoretically expected, the power law part of spectra is almost independent of the soft 
photons distribution. In other words, it is independent of the illumination of the bulk 
motion site by the accretion disk flow. 


In Figure 5 we give the result of our computations for kT^. = 20 keV, and m = 4 which 
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reproduces the main features of the spectra observed in the high state of BHS: a spectral 
index of 1.5 and an high energy cutoff around 500 keV. The effect of the thermal motion 
for temperatures higher than 10 keV is there clearly seen. The spectrum is getting harder, 
i.e. the index decreases and the high energy cutoff increases. 

Furthermore, the spectra obtained at high kT^, and m greater than 2 have a spectral 
index around 0.8 which is what is observed from BHS in the hard state, where the thermal 
Comptonization is thought to dominate (ETC). 

This is also shown in Figure 6 where we present the results of the computation for 
kTf. = 50 keV, and m = 4, compared to the analytical solution of Hua & Titarchuk 1995 
(Eq. 6) derived for the pure thermal Comptonization case. The spectrum is practically 
identical to the standard thermal Comptonization spectrum. In this case the effect of the 
bulk motion can be seen only at high energies where there is an excess in the converging 
inflow spectrum due to coupling of the thermal and bulk motion velocities. It is worth 
noting here that this excess is really detected in the observation of the hard state of Cyg 
X-1 source (e.g. Ling et al. 1997, see also for details of the model, Skibo & Dermer 1995). 

Editor, please put Figure 6 here. 

So, the change of spectral shapes from the soft X-ray state to the hard X-ray state 
is clearly to be related to the temperature of the bulk flow. In other words the effect of 
the bulk Comptonization compared to the thermal one is getting stronger when the plasma 
temperature drops below 10 keV. Also, as it could be seen by comparing the results in the 
“Flat” and general relativistic case, the spectral index of the Comptonized spectrum is very 
dependent of space curvature effects, and then its observation should be a very powerful 
tool to determine the nature of the central object. 
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5. Spectral State Transition and Temperature of the Converging Inflow 

Atmosphere 

To complete the discussion of the Monte Carlo simulations and their application to the 
observed spectra of BHC systems we present some details of the thermal balance treatment 
of the bulk motion atmosphere (cf. Chakrabarti & Titarchuk 1995, Titarchuk, Lapidus & 
Muslimov 1998). 

The total count rate from the source increases with an increasing accretion rate. The 
X-ray spectrum becomes softer, due to the increase of the supply of soft photons from the 
disk illuminating the advection dominated region (shocks, centrifugal barrier) where the 
energy release is due to geometrical compression of the non-Keplerian accreting matter. 
This region, the so called Compton cloud, is the source of hard X-ray radiations which are 
produced as a result of the upscattering of the soft photons emitted from the disk. In the 
flnal stage of the accretion (this is within a few Schwarzschild radii) the matter goes towards 
the black hole almost free falling as the radiation pressure cannot stop matter falling in. 

The Compton cloud (CC) region can be treated as a potential wall at which the 
accreting matter releases its gravitational energy. This occurs in an optically thin region 
where the column density is of order of a few grams, or where the Thomson optical thickness 
To ~ a few. The amount of energy released per second is a fraction of the Eddington 
luminosity since the CC region is located in the very vicinity of a central object (~ 3-6 
rg). The heating of a gas due to the gravitational energy release should be balanced by the 
photon emission. For high gas temperatures, Comptonization is the main cooling channel, 
and the heating of electrons is due to their Coulomb collisions with protons. Under such 
physical conditions the energy balance can be written as (see e.g. Zel’dovich & Shakura 
1969, hereafter ZS69, equation [1.3]) 


F/to = Co-^{a)e{T)Te/f{Te). 


( 3 ) 
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Here r is the current Thomson optical depth in the emission region (e.g. in a slab), a is the 
energy spectral index for a power-law component of the Comptonization spectrum, e{T) is 
a distribution function for the radiative energy density, f(Te) = 1 -|- 2.5(kTe/meC^), Tg is 
the plasma temperature in K, Co = 20.2 cm s“^ is a dimensional constant, and, hnally, 
99 ( 0 ;) = 0.75q;(1 -|- q;/ 3) if (p{a) < 1, otherwise = 1. The latter formula is obtained by 
using the relationship between the zero- and first-order moments (with respect to energy) 
of the Comptonized radiation field (Sunyaev & Titarchuk 1985, §7.3, equation [30]). The 
distribution of the radiative energy density in the emission region e(r) can be obtained from 
the solution of the diffusion equation (cf. ZS69, equation [1.4]), 

ld?e _ _ (F/c) 

3 To 

subject to the two appropriate boundary conditions. 

The hrst boundary condition must imply that there is no scattered radiation from the outer 
side of the emission region, i.e. 

H S 

--= 0 for r = 0. (5) 

dr 2 

In our case this condition holds at the inner surface of a slab facing a central object (black 
hole). In the neutron star case one can expect some additional soft flux from a NS resulting 
in an additional illumination of the inner surface of a slab. It is easy to generalize the 
following analysis for this particular case and we offer the reader to make this exercise. 


The second boundary condition requires that at the outer surface of a slab the incoming 
flux should be equal to the external flux i.e. 


1 de 
3 dr 


c 


at T = Tq. 


( 6 ) 


The solution of equations (4)-(6) provides us with the distribution function for the 
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energy density 


£{r) = + 3to[t/to - 0.5(r/ro)^F/(F + Fd)]}. 


Thus, from equations (3) and (7) we get 


ip{a)T^To 

fiTe) 


^ 0.75 ■ 10^ 


F 

F + Fd 


K. 


( 7 ) 

( 8 ) 


When Fd F the spectral index a varies very little since a is a function of Tero/f(Te) (see 
e.g. Titarchuk & Lyubarskij 1995 for the thermal Comptonization case). Thus, as long 
as the external flux (due to the photons from the disk) is much smaller than the internal 
energy release (per cm^ per second) in the CC region, the spectral index is insensitive to 
the accretion rate in the disk. The values of parameters consistent with equation (8) are: 
To ^ 5, Tg ^ 2 X 10® K, and a ^ 1, which are characteristic of a hard state for the galactic 
BH and NS systems. 


When the external flux Fd becomes comparable to the internal energy release, F, the 
cooling becomes more efficient not only due to Comptonization, but also due to the free-free 
cooling, and therefore the electron temperature unavoidably decreases, Tg 10®-F/(F-|-Fd) 
K (see CT95 for the numerical calculations of spectral indices and temperature). 

This illumination effect can explain the hard-soft transition when the temperature of 
the Compton cloud drops substantially with the increase of the soft photon flux from the 
disk. As it is seen from above (Eq. 8) the temperature drops from 50-60 keV in the hard 
state (i.e. when Fd -C F) to 5-10 keV in soft state (when Fd and F are comparable). For 
this particular range of temperatures the Monte Carlo simulations have been done in this 


paper. 



20 


6. The relevance of Bulk motion Comptonization model to the recent high 

energy observations 

Shrader and Titarchuk, (1998), have already shown that the spectral shape of Compton 
scattered photons in converging flow is very close to the observed spectra. However, it is 
still necessary to show how many photons escape from the central region, when a given 
number of photons are injected in the Monte-Carlo simulations. For example in the case 
of kTe = 20 keV, and m = 4 (see Fig.5), among 10^ generated photons there are 15,372 
photons, or 15.4 %, escaping without interacting, 45,557 photons, or 45.6 %, going through 
the black hole horizon, 13,393 photons, or 13.4 %, scattered in the accretion disk (but 
not in the converging Compton cloud), the inner radius temperature of which = 0.5 keV, 
and 25,678 photons, or 25.7 %, scattered in the Compton cloud (some of these photons 
being also scattered in the disk). Thus almost the same photon flux (for our particular 
geometry, see Fig. 2) 28.8 %, (15.4 % -|- 13.4 %) goes directly to observer and 25.6 % after 
scattering in the converging Compton cloud. But a smaller flux of the detected photons 
emerges after multiple scattering and they form the hard tail of the resulting spectrum 
(Fig. 5). We remind the reader (see for details ST80, T94, and TMK97) that the photons 
undergoing multiple scatterings produce the specihc space distribution (in accordance to 
the hrst space eigenfunction, see TZ98, Fig. 3) and that fraction of photons undergoing 
multiple scattering in the plasma cloud, fms is related to the expansion coefficient of the 
space source distribution over the hrst space eigenfunction. In the simplest case, the 
uniform source distribution fms is approximately 0.8. But in our case when the source 
photon distribution is produced by an external illumination of the converging inhow (Cl) 
region by the soft disk photons and the optical depth of the converging region is of order 
of 1 (r ~ m[1.5“^/^ — 3“^/^] ~ 1), this fraction is slightly less. Here for the optical depth 
estimate we use formula (2) with an assumption that rout = and = l.br^ being as a 
photon bending radius (TZ98). 



21 


Thus the fraction of photons undergoing multiple scattering and detected by the 
observer with respect those photons undergoing a few scattering and detected or comes 
directly to the observer from the disk is estimated as follows 




0.26-0.8 
0.26-0.2 + 0.29 


0 . 6 . 


(9) 


The exact number of these photons can be determined through the htting of the 
simulated resulting spectrum to the test spectrum, composed of two components: a 
Planckian spectrum with a hxed temperature of the soft photons and convolution of some 
fraction / of this component, with the Cl upscattering Green’s function. It is easy to 
prove (using the property of the Green’s function to conserve the number of the scattering 
photons) that / is exactly the fraction of photons (with respect to the thermal component) 
multiply scattered in the converging inflow. As shown in §3, / = 0.48 (for our case) and the 
flux of the escaping hard photons is large enough to explain the observations. We remind 
the reader that the empirically determined fraction are 0.32 and 0.72 for GRO J1655-40 
and GRS 1915+105 respectively (ShT98) and thus Comptonization by the converging inflow 
is consistent with the observations. 

It is worth noting that the number of the photons forming the very hard tail (for 
energies E 3> hTfl) 

Nhard{E)= r I(Eo)^ ( 10 ) 

Je Eq 

is much less than the number of photons in the thermal component N^h, i.e. 

< f{E/2kTfl-‘^ < 1 . ( 11 ) 

Nth 

For example, this ratio is about 8 - 10“^ for E = 10 keV in the cases of GRO J1655-40 
and GRS 1915+105. In other words only 8% of the photons seen as coming from the disk 
are effectively scattered off electrons within a few Schwarzschild radii (2 — 3 Vg) and form 
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the extended power law tail of the spectrum. Inequality (11) can be readily proved using a 
formula for the hard tail energy flux I{E) (Eq. 10) presented in the convolution form (see 
ShT98, Eq. 2). 

In the end of this section we discuss a issue related with the illumination of the Cl 
atmosphere by the disk photons. For the beginning we assume that the disk emitting 
isotropically is extended from the inner radius to the outer radius rout and the outer 
radius of the Cl atmosphere is r^p. And further more we assume that Vgp < Tin- 

In fact, the disk flux intercepted by the the Cl atmosphere with radius Vgp (see Fig. 2) 
is calculated through the integral 

front fTt/2 fl 0 1/0 

Pint = TT / rdr / cos ifdif / (1 — /i ) ^^dfj, (12) 

Jrin 2-7r/2 J Pt 

where /i* = — {rsp/ry. In the case when the disk illuminates the slab (with the half-width 

Ar = Tin = Sts and the half-height H) situated in the center /x* = {r—Vin)/\j— VinY + H'^. 

The two internal integrals of Eq. (12) are calculated analytically. The first one is equal 
to 2 and the second one is 6^* — 0.5 sin26**, where 6** = arccos/x*. Even the triple integral 
of Eq. (12) can be estimated analytically with accuracy better than 10% for the spherical 
case with the aforementioned assumption that Vgp < Tin- 

Pint = Y^lpir--^ - r;J)). (13) 

In the case when the part of the disk is covered by the Cl atmosphere the flux. Pint 
is calculated by integration from r-gp to Tout followed by adding 7r^(r^p — ffn)-, the disk flux 
emitted inside the atmosphere. 

The fraction of the disk photons illuminating the Cl atmosphere is 

PintPsc 

T^'^irLt-rfn)' 


fill 


(14) 
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Where Psc is a probability of scattering of the disk photons in the Cl atmosphere which is 
of order 1 — e~'^. 


The factor /, used in the spectral htting can be estimated from here 


( 15 ) 


{I - fms)fillA + p{l - fill) 

where /i is the cosine of inclination angle of the disk and A is a spherical albedo of the Cl 
atmosphere. 


The quantity A can be estimated as follows (e.g. Sobolev 1975) 


A 


3t/4 + 1 ■ 


( 16 ) 


This formula is obtained in the Eddington approximation with an assumption of isotropic 
scattering at any event and the pure absorptive inner boundary. Despite of the Eddington 
approximation the albedo formula (16) provides an accuracy of order 10% for any optical 
depth. Eor optical depth 1, A = 0.42 which is slightly higher than value obtained in our 
Monte Carlo simulations, Amc = 25.7/(25.7 + 45.6) = 0.36 (see above) because of the 
photon bending. 


Now we present the numerical estimates of the factor / for two different cases. The 
Erst one related with our Monte Carlo simulations when ~ 1 — e~^ = 0.63, fms = 0.8, 
A = 0.42 and we assume that /i ~ 1. Then we get that / ~ 0.5 which is very close to that 
we get from the spectral htting of the Monte Carlo simulated spectrum. It is worth noting 
here that in the limit of the high optical depth (or mass accretion rate) / ~ /ms/(l — fms) 
because psc, fui, and A converge to one. Thus in this particular case the factor / can be of 
order of one or more. 


For another case we assume that the soft radiation comes from the extended disk with 
the inner radius = Svg and Tout = which is the size of the disk emission region 
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related with the color temperature of 0.85 keV (for details of this model see Borozdin et 
ah 1998) and the inclination direction cosine /i = 0.5. Also we use the above values of 
Psc ~ 0.63, fms ~ 0.8, A ~ 0.42. Using Eqs (12), (14-15) we obtain that / = 0.68 ■ 10“^ if 
we consider the illnmination of the sphere sitnated in the center with Vgp = Sr*. The factor 
/ is 1.7 ■ 10“^ in the case of the illnmination of the slab with the half-height, H = Sr^. 

To conclnde this section we emphasize once again that the Comptonizing region must 
overlap with the inner region of the disk, in order to get the high fraction of the seed 
photons from the disk (see ShT98). Bnt from the other hand the high fraction, / (of order 
10 % and higher) can be also obtained if the innermost region of the disk is pnffed np 
to heights of order of 3 r* and the disk is almost seen at edge on. More details of this 
discussion and comparison with the RXTE observations would be presented in onr incoming 
paper (Borozdin et al. 1998). 


7. Conclusions 

We show that spectra in the soft and hard states are formed by npscattering of the 
disk soft photons by the matter (electrons) moving in the strong gravitational held. 

The presence of the event horizon as well as the behavionr of the null geodesics in its 
vicinity largely determine the dependence of the spectral index on the how parameters. The 
efficiency of gravitational to radiative energy conversion is of order 5% or less in the disk 
aronnd black hole (e.g. SS73; Narayan, Barret & McClintock 1997). Thus it is natural to 
expect that radiation pressnre cannot hold back the matter which is almost in free-fall. Bnt 
the shape of the spectrnm (the spectral index and the position of the high energy cntoh) 
depends on the temperatnre of the bnlk motion atmosphere, which is determined by the 
efficiency of the cooling of the atmosphere by the flnx coming from the disk (see §5 here 
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and CT95; Titarchuk, Lapidus & Muslimov 1998). In the hard state the most gravitational 
energy is released out of the disk (see possible scenarios in CT95) and the bulk motion zone 
is getting very hot [because there is no enough cooling agents (soft disk photons) to cool 
down the environment]. The plasma temperature is of the order of 50 keV. 

In this case, the contribution of the outer part of the bulk motion atmosphere in the 
resulting spectrum is signihcant because of the large volume of this region and of the high 
thermal velocities. In fact, the hard photons are produced mainly by effective scatterings 
of the disk soft photons on electrons at angles near tt (Hua & Titarchuk 1995, Fig. 2). 
But in the outer layer the gravitational bending and redshift is less pronounced and thus 
the effect of General Relativity on the spectral index is smaller for larger values of electron 
temperature as shown in table 2. Indeed, the spectrum is practically identical to the 
standard thermal Comptonization spectrum (see Fig. 6). 

When the mass accretion rate in the disk increases, the plasma temperature goes down 
to such a level (of order 10 keV or less) that only the free-falling electrons transfer their 
momentum to the soft-photons producing the power-law component extending to energies 
comparable to the kinetic energy of electrons in the converging inflow, i.e. of order rrieC^- In 
this case, the spectrum observed at inhnity consists of a soft component produced by those 
input photons that escape after a few scattering without any signihcant energy change 
and of an hard component (described by a power law at the energies much higher than 
the characteristic energy of the disk soft photons) produced by the photons that undergo 
signihcant upscattering. 

The luminosity of the power-law component is relatively small compared to that 
of the soft component. The power law is quite steep due to the strong photon bending 
and the gravitational redshift near the black hole (within 1.5-2 Schwarzschild radii). We 
demonstrate through our Monte Carlo simulations how many photons escape from the 
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central region, when a given number of photons are injected in the Monte-Carlo simulations 
and that fraction is perfectly consistent with the recent high energy observations (ShT98). 

We also show that our Monte-Carlo spectra are htted quite well by the the analytical 
bulk motion Comptonization model (BMC) which allows one to use this model for the 
htting of the real observational data. The analytical model (BMC) contains only three 
parameters. Two of them are the color temperature of the seed photons Tc and the 
relative weight of the soft component in the entire spectrum, (1 — /) which are the main 
characteristics of the soft component. Thus it is getting possible to determine the color 
temperature dependence on the energy flux in the soft energy band. The disk origin of the 
soft seed photons can be conhrmed or refuted by this dependence (see details in Borozdin 
et ah 1998). 

We demonstrate using our Monte-Carlo simulation that the main features of the 
simulated spectra ( a spectral index of 1.5 and an high energy cutoff around 500 keV) 
are consistent with the observed X-ray spectra in the high state of BBS (e.g. Grove et 
ah 1998). The predicted high energy cutoff (assuming a relativistic Maxwellian electron 
distribution) is in the energy range between 200 and 400 keV (see Figs 4-5). It is slightly 
less than that observed in the superluminal source GRO J1655 by OSSE (Grove et ah 
1998, table 2 and Fig. 2). It is worth noting that the break energy is quite sensitive to the 
electron distribution. In future work, we will show that the observed spectrum (including 
the position of the high energy break) can indeed be accounted for by the converging 
model quantitatively using other electron distributions. In fact, we plan to investigate the 
possibility of non-thermal tails to the electron distribution in accreting systems and the 
expected radiation signatures in the hard X-ray band. 

Thus our results of the Monte Carlo simulations strongly support the idea that the bulk 
motion Comptonization might he responsible for the extended power-law speetra seen in the 
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black-hole X-ray sources in their soft state. And the hard-soft states transition is regulated 
by the plasma temperature of the converging inflow into a black hole. 
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Fig. 1.— The left-hand panel is an artists conception of an accretion disk surrounding a 
black hole. The approximate sites from which the soft-seed radiation is produced, the more 
general soft radiation which is directly viewed by the observer (thick solid line), and the 
hard-radiation emission ((thin solid line) are indicated. Also, the dotted lines indicate the 
scale and location of the bulk-motion scattering site. The schematic diagram in right-hand 
panel, depicts a typical sequence of scatterings within the bulk-motion inflow region. Some 
details are given in the text. 

Fig. 2.— Scheme of the geometrical setup used in the Monte-Carlo simulations. We vary 
the accretion sphere radius rout from 3 to 20 in our simulations (see text). 

Fig. 3.— The Monte Carlo photon spectrum (see below Fig. 5) is htted to the test 
(analytical) spectrum, composed of two components: a Planckian spectrum with a hxed 
temperature of the soft photons and convolution of some fraction of this component with 
the Cl upscattering Green’s function. 

Fig. 4.— Emergent Comptonized photon spectrum in the general relativistic case: in this 
simulation, the cloud has a temperature of 5 keV and m is equal to 2. The best power law 
fit also shown has a a photon number index of 2.9 (or the spectral index a = 1.9). 

Fig. 5.— Emergent Comptonized photon spectrum in the general relativistic case: the 
physical conditions of the cloud are now: kT^ = 20 keV, m = 4. The best power law £t also 
shown has a photon number index of 2.5 (or a spectral index a = 1.5). 

Eig. 6.— Emergent Comptonized photon spectrum in the general relativistic case with this 
time kTe = 50 keV, and m = 2. The high energy part of the spectrum is compared to the 
thermal Comptonization spectrum [Hua & Titarchuk (1995), Eq. 6]. 
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